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(Dated: June 2, 2008) 

Abstract 

Using 281 pb~^ of data collected with the CLEO-c detector at the ^(3770) resonance, wc have 
studied Cabibbo-suppressed decays of D mesons to final states with two kaons. We present results 
for the absolute branching fractions of the modes K~^K~ , — > K^Kg, and K~^Kg. 

We measure B{D^ K+K') = (4.08 ± 0.08 ± 0.09) x 10-^ B{D^ K^^^) = (1-46 ± 0.32 ± 
0.09) X 10"^, and B{D+ K+Kfj = (3.14 ± 0.09 ± 0.08) x lO'^. We also determine the ratio 
^(L>° K+K-)/B{D^ TT+vr-) = 2.89±0.05±0.06. For each measurement, the first uncertainty 
is statistical and the second uncertainty is systematic. 
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Detailed studies of the decay of charmed mesons provide an opportunity to probe the 
interplay between the strong and weak force. One such exarnple is the study of the contribu- 
tions of l^-exchange diagrams and final state interactions [jj. In particular, the branching 
ratio B{D^ K^K^)/B{D^ tt+tt^) exposes large SU(3) flavor symmetry breaking ef- 
fects. To first order one would expect a branching ratio close to unity, but the most recent 
studies show the ratio to be 2.760 ± 0.040 ± 0.034 \2\. The direct weak decay K^K^ 
is expected to occur primarily via two W^-exchange diagrams, which in a four-quark model 
under SU(3) is expected to be zero due to destructive interference 0. Hence, substantial 
contributions to the rate for K^K^ are expected from rescattering in such channels 

as L>° ^ K*+K*~, tt+tt", and D° ^ p+p~ 0, H, Q- Theoretical estimates of SU(3)- 

breaking effects give branching fractions ranging from a few parts in 10^ [Hj up to 0.3% 
Here we study the decay mode —>■ K^Kg, whose branching fraction should be half that 
of — > K^K^. Note that the decay mode K^K^ is forbidden by CP conservation. 

Understanding SU(3) flavor symmetry breaking is important in attempting to quantify Stan- 
dard Model contributions to — mixing 0. The decay mode D+ K+K^ is useful 
for the estimation of SU (3)- violating effects in the D meson system [8|. Thus, a precise 
measurement of the absolute branching fraction of D — > KK modes will test our current 
understanding of the Standard Model. 

We use a 281 pb~^ data sample of e^e~ collisions at a center-of-mass energy of approx- 
imately 3770 MeV produced at the Cornell Electron Storage Ring and recorded with the 
CLEO-c detector. The CLEO-c detector is a general purpose solenoidal detector that in- 
cludes a tracking system for measuring momenta and specific ionization {dE/dx) of charged 
particles, a ring imaging Cherenkov detector for particle identification, and a Csl calorime- 
ter for detection of photons. The details of the CLEO-c detector have been described 
elsewhere |i, Eol, Ell • 

We are interested in events of the type e+e~ — > '?/'(3770) DD. The D mesons have 
energy equal to the beam energy since the ?/'(3770) resonance is below the kinematic thresh- 
old for DDtt production. After selecting candidate particles that may be the daughters 
of a D meson decay, we define the following two variables, Ai? = Ylii ~ -^'beam and 
Mbc = \/ -^beam ~ I Si PiP' wherc Ei and Pj are the energies and momenta respectively of 
the D meson daughters and £'beam is the beam energy. If the event is reconstructed cor- 
rectly, AE should be consistent with zero and the beam constrained mass Mbc should be 
consistent with the mass of the D meson. 

In selecting events, charged tracks are required to satisfy criteria based on the track 
quality [13]. They must also be consistent with coming from the interaction point in three 
dimensions. Pions and kaons are identified by consistency with the expected dE/dx and, 
when available, ring imaging Cherenkov information. We also require that |Aii^| be less than 
20 MeV. Other selection criteria are mode dependent as described below. Charge- conjugate 
modes are imphed throughout. 

In selecting events for the mode K^K~ , we veto events that contain leptons, e.g. 

cosmic rays and Bhabha events, and furthermore require cos 6* > —0.9 for K" and cos 6* < 0.9 
for K~^, where 6 is the angle between the e+ direction and the kaon momentum direction. 
The cos^ criterion is useful to further suppress radiative Bhabha background events that 
have the same energy loss at about 800 MeV as the kaons in D° ^ K+R-. 

In selecting events for the mode D'^ — > K^K^^ we require the candidate invariant 
mass to be within 12 MeV/c^, or roughly 4.5cr, of the known invariant mass 13j and 
require the candidate to have a flight distance significance {L/ai) greater than two in 
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FIG. 1: Mbc plots for data, left to right: L>° K+R- and 
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FIG. 2: K^K^ plots for data, left to right: Mbc fits for the signal region for K^s^^s^ 

sideband region 1 for — > K^K^ and sideband region 2 for —> K^Kg. The sideband regions 
are shown in Fig. [3l 

order to separate true Kg particles from pairs of pions produced at the interaction point. 

The mode —>■ K^K^ has significant backgrounds, mostly from the mode 
Kg7T~^7r~ where the vr"'" and the 7r~ have an invariant mass consistent with the Kg mass. 
Therefore, similar to our analysis of the mode K~^Kg, we require the flight distance 

significance for each of the K^ candidates to be greater than two. Furthermore, we perform 
a background subtraction where we take the region |M^+^- — M^o| < 7.5 MeV/c^ to be our 
signal region and take the region 20 > \M^+^- — M^o| > 9 MeV/c^ to be our background re- 
gion, where M^^+j^- is the reconstructed Kg candidate invariant mass and Mj^o is the known 
Kg invariant mass [is']. This will be described in detail later. 

The branching fractions are calculated using B = Y /2eNj^[), where Y is the measured 
signal yield, e is the efficiency obtained from Monte Carlo simulations, and Nj^q is the 
produced number of DD pairs (charged or neutral depending on the decay mode) in our 
data sample. We use A''^,^ from Ref. |l2l |. 

To determine the measured yields, we perform an unbinned maximum likelihood fit to a 
signal shape and one background component to extract signal yields from the Mbc distribu- 
tions. The signal shape includes the effects of beam energy spread, initial state radiation, the 
line shape of the V^(3770), and reconstruction resolution 12]. The background is described 



by an ARGUS function [ij], which models combinatorial contributions. We have validated 
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TABLE I: Efficiencies and yields with their statistical uncertainties. The yield for D KgKg is 
the yield after background subtraction. 





Mode 


Efficiency (%) 


Yield 




^ K+K- 


56.47 ±0.33 


4746 ± 74 






22.75 ±0.22 


68 ± 15 


D+ 




38.29 ±0.28 


1971 ± 51 
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FIG. 3: (M^+^- - Mj^o) vs (M^+^- - M^o) for D° ^ ^5^5 candidates. The region labeled "S 
is the signal region, the regions labeled "1" are the sideband 1 region, and the regions labeled "2" 
are the sideband 2 region. 

our analysis procedure by studying simulated DD events in a sample 40 times the size of our 
data set and we reproduce the input branching fractions to within statistical uncertainties. 
Our fits to the data are shown in Figs. [1] and [21 and our efficiencies, along with the fitted 
yields, are given in Table [B 

The scatter plot for the two candidate invariant masses in data in the region |Mn-+^- — 
Mf^o\ < 20 MeV/c^ is shown in Fig. [3l where M^+tt- is the reconstructed candidate 



invariant mass and M^o is the known Kg invariant mass 13|. The order of the two Kg 
candidates is essentially random as it is determined by the order the tracks were found 
in the event. We expect the bulk of the signal to be in the central (signal) region which 
we define to be |M^+^- — M^o | < 7.5 MeV/c^ for the two Kg candidates. In order to 
account for K^-k^ti^ and other background events entering into the signal region of 

KgKg, we perform a background subtraction. We then take the four sideband regions 
where one Kg candidate has 20 > \M^+^- — Mp^o \ > 9 MeV/c^ and the other Kg candidate 
has IMtt+tt- ~ Mj^o \ < 7.5 MeV/c^. We call this region sideband 1. Finally we take the four 
corner regions where both Kg candidates satisfy 20 > {M^^+t^- — Mj^o \ > 9 MeV/c^. We call 
this region sideband 2. The number of candidate events in each region is determined by a 
fit to the Mbc distribution of the D meson candidates and are found to be 5* = 95 ± 13, 
Si = 40 ± 11, and 6*2 = 2 ± 3 for the signal region, sideband 1, and sideband 2 respectively. 
We let ri denote the ratio of the area in the signal region to the sideband region 1 and r2 



5 



TABLE II: Measured branching fractions and Particle Data Group (PDG) values 13|]; the first 
uncertainty is statistical and the second uncertainty is systematic. 





Mode Our Result 




PDG 2007 


B{D+ 


K+K~) (4.08 ± 0.08 ± 0.09) x 
K'^s^s) (1-46 ± 0.32 ± 0.09) x 
K+K°) (3.14 ± 0.09 ± 0.08) x 


10-3 
10-4 
10-3 


(3.85 ±0.09) X 10-3 

(3.6 ±0.7) X lO"'' 
(2.95 ±0.19) X 10-3 



denote the ratio of the area in the signal region to the sideband region 2. Then the total 
yield Y is given hj Y = S — 2riSi +r2S2- Note that Si is the measured yield in two bands, 
so after rescaling 5*1 with ri, a factor of two is needed in the second term on the right hand 
side to account for the fact that two bands are overlapping in the signal region. The third 
term is added since the second term oversubtracts the uniform background, but we find that 
5*2 is much smaller than S and 5*1. We have ri = 225/660 and r2 = 225/484. 

We have studied various systematic uncertainties in our analysis. These include particle 
identification efficiency, tracking efficiency, uncertainties in the fitting procedure, and the 
uncertainties due to the AE selection criterion. Using techniques outlined in Ref. fl2| . 
the systematic error in the absolute tracking efficiency for pions (kaons) was determined 
to be 0.3% (0.7%). We have also determined that our Monte Carlo simulation accurately 
models the identification efficiency of pions, kaons, and Kg to within 0.3%, 0.2%, and 1.9% 
respectively. To determine the systematic errors due to the fitting procedure, we varied the 
background shape and the momentum resolution in the fit and compared the variations in 
the fit yield. These have been determined to be 0.8%, 0.6%, and 4.1% for each of the modes 
— > K~^K~, — > K~^Kg, and — > K^Kg respectively. To determine the systematic 
error due to the AE criterion, we varied the AE criterion and studied its effect on the 
yield. This has a negligible effect for the modes K~^K~ , —>■ K~^Kg, and for the 

mode —>■ K^Kg we assign a systematic error of 2.0%. Finally, using a detailed Monte 
Carlo simulation of the mode —>■ Kgir'^Tr' , we estimated the effect of non-linearity in the 
background shape for the background subtraction in the mode —>■ KgKg and found it to 
be negligible. 

The measured branching fractions are summarized in Table [Tll We use the DD yields 
Njjono = (1.031 ± 0.008 ± 0.013) x 10^ and Nd+d~ = (0.819 ± 0.008 ± 0.010) x 10^ as 
measured in Ref. [12] to determine our branching fractions. We also determined the ratio 
of branching fractions B{D^ — >■ K^K~) to B{D^ tt^tt~), where the latter branching 
fraction has previously been measured using the same data sample [ll], to be i3(D° — >■ 

K") / B{D^ vr+TT-) = 2.89 ± 0.05 ± 0.06, consistent with the previous measured value 
of B{D^ K+K-)/B{D° vr+vr-) = 2.760 ± 0.040 ± 0.034 0. The modes D° ^ K+R- 
and — > K~^Kg are consistent with the PDG average. Finally, for the — > KgKg 
mode, our errors are significantly smaller than previous measurements, while our branching 
fraction is about a factor of 2.4 (2.7 standard deviations) smaller (l3| . 
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